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A model Vlasov-Poisson system is simulated close the point of marginal stability, thus assuming 
only the wave-particle resonant interactions are responsible for saturation, and shown to obey the 
power-law scaling of a second-order phase transition. The set of critical exponents analogous to 
those of the Ising universality class is calculated and shown to obey the Widom and Rushbrooke 
scaling and Josephson's hyperscaling relations at the formal dimensionality d — 5 below the critical 
point at nonzero order parameter. However, the two-point correlation function does not correspond 
to the propagator of Euclidean quantum field theory, which is the Gaussian model for the Ising 
universality class. Instead it corresponds to the propagator for the fermionic vector field and to 
the upper critical dimensionality d c — 2. This suggests criticality of collisionless Vlasov-Poisson 
systems as representative of the universality class of critical phenomena of a fermionic quantum 
field description. 
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I. INTRODUCTION 

The remarkable property of critical phenomena is the 
universal scaling appearing in vast variety of systems; 
e.g., magnets and gases follow simple power laws for 
the order parameter, specific heat capacity, susceptibil- 
ity, compressibility, etc. [l]j. In thermodynamic systems, 
phase transitions take place at a critical temperature T cr 
when the coefficients that characterize the linear response 
of the system to external perturbations diverge. Then 
long-range order appears, causing a transition to a new 
phase due to collective behavior of the entire system [!| • 

The condition for nonlinear saturation in the test case 
of the bump-on-tail instability in plasmas [j| is 

uj b « 3.2 7L (1) 

0, where 7l is the linear growth rate of a weakly 
unstable Langmuir wave according to the Landau the- 
ory 6], u>b — [eEk/m) 1 / 2 is the frequency of oscillations 
for particles trapped by the wave. These trapped par- 
ticles generate a long-range order of the wavelength fc, 
the saturated amplitude E can be considered as the or- 
der parameter, and the condition of saturation can be 
rewritten as a power law, typical for the second-order 
phase transitions, E ~ 7^. 

However, unlike thermodynamics this scaling contains 
the nonthermal control parameter 7 £, which is deter- 
mined by the slope of the distribution function dfo / dv 
at the phase speed v r — 0J pe /k of the perturbation near 
the electron plasma frequency u> pe . [For thermodynamic 
systems like magnets for which the magnetization M be- 
low the Curie point T cr is the order parameter, the scal- 
ing is M ex. e@ , where e = (T — T cr )/T cr at T < T cr .] 
Another difference is the critical exponent itself - rela- 
tion Q predicts the very unusual exponent 8 = 2; in 
contrast, the hydrodynamic Hopf bifurcation |7fl, also de- 
scribed by the same scaling between the saturated am- 
plitude and the growth rate, has the mean-field critical 
exponent 8 — 1/2. 



An analysis [8( assuming thermalization in a Vlasov- 
Poisson plasma or gravitating system leads to a critical 
exponent 8 < 1, and the exponent 8—1/2 has also been 
hypothesized for the bump-on-tail instability in [9j] . How- 
ever, detailed center-manifold analysis, which establishes 
the normal form for a weakly unstable perturbation in 
one-component collisionless Vlasov-Poisson system, con- 
firms 8 = 2 llCjl . The exponent 8 = 2 is also confirmed 
numerically |lll IT^ . 

The striking discrepancy between these exponents can 
be better understood if we consider the structure of the 
phase space corresponding to these cases. The exponents 
8 =1/2 and (3 < 1 Q correspond either to saturation 
of a strongly dissipative instability or to a thermalized 
system. In both these cases the distribution function can 
be factorized as f(q,p) = y(q)g(p), where g(p) can be as- 
sumed to be Gaussian, and the system is described by its 
momenta. The exponent 8 = 2 corresponds to saturation 
due to nonlinear wave-particle interactions in a weakly 
unstable collisionless system, where correlations between 
coordinates and impulses are not destroyed by dissipa- 
tive processes, so the description cannot be reduced to 
moments of the distribution. 

More formally, a dissipative and/or thermalized sys- 
tem is represented by a discrete set of momenta of the 
distribution function f(x, v, i), which depend only on the 
coordinate x, but not on the velocity v, M = {p,v,T}, 
where p, v, and T are the local density, the velocity, and 
the temperature, respectively. The evolution is a flow 
g t which maps M onto itself, g t : M — > M. In fact, in 
a neighborhood of the threshold 7 l — the evolution 
gt'- M — > M can be reduced to a normal form, which 
maps only the order parameter, g' t : n —> n, where n = M° 
(or C , where C is the set of complex numbers), and 
therefore the evolution is a trajectory n = n(t) or in 
other words the set Y = M + x R°. The phase space 
of a one-dimensional collisionless system is a continuous 
set JI = 1x1 and evolution can be represented as the 
flow Wt'. H — > H. (For periodic boundary conditions the 
phase space is isomorphic to a cylinder C = Txl, where 



2 



T is isomorphic to a circle.) The sets n and H (or C) have 
different dimensionality, and therefore renormalization of 
collisionless system in a vicinity of threshold - i.e., trans- 
formation of the set H and the mapping wt - involves one 
or more additional dimension. Further it is shown below, 
that scaling transformations close to the threshold are in- 
terrelated with the additional velocity coordinate, which 
disappears in hydrodynamic description because of inte- 
gration of the distribution function f(x, v, t) over v. 

From the theory of critical phenomena it is known that 
dimensionality d is inseparable part of the threshold de- 
scription - along with the critical exponents (e.g., [l3|p. 
Besides f3, other critical exponents: a, 7, S, v, and 77, 
describe the following scalings of the Ising universality 
class: (i) the specific heat capacity scales as 



(ii) the susceptibility as 



X 
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(iii) the response M at e = as 
M (x B 1/s 



cx e 



(2) 



(3) 



(4) 



(iv) the correlation length as 

£~N-"; (5) 

and (v) the two-point correlation function as 



e -r/€ 

G ( r ) ~ r d-2+r. 



(6) 



These exponents are not independent, but are interre- 
lated via scaling laws, e.g., the Widom equality 



7 = /?(<5-l) 



(7) 



|14| . These scaling laws also include hyperscaling laws 
such as Josephson's law, 



vd 



a , 



(8) 



|l5l | which involves the dimensionality d along with the 
exponents. 

For thermodynamics the mean-field exponents are of 
the Landau-Weiss set, a = 0, /3 = 1/2, 7 = 1, 8 = 3, 
v = 1/2, and rj — 0, and the scaling laws hold at the for- 
mal dimensionality d = 4. However, the possibilities of 
critical phenomena are not exhausted by the Ising uni- 
versality class - the percolation critical exponents [T^ |. 
which describe another vast class of critical phenomena, 
are different from those in thermodynamics and scaling 
laws hold at a different dimensionality. In particular, for 
the Bethe lattice (or Cayley tree) [T3 Josephson's law 
holds at dimensionality d = 6. Despite the description 
being the same, this difference separates the cases into 



different universality classes with different upper critical 
dimensions: d c — 4 for the Ising universality class |l8j 
and d = 6 for percolation. 

For a collisionless gravitating system, where the sat- 
uration mechanism is the same as for the bump-on- 
tail instability in plasmas, the critical exponent (3 = 
1.907±0.006, and the critical exponents 7 = 1.075±0.05, 
5 = 1.544± 0.002 can be determined analogously to ther- 
modynamics and calculated from the response to an ex- 
ternal pump [l9| . These exponents are very different 
from the thermodynamic set, but nevertheless satisfy the 
Widom equality, thus suggesting the validity of scaling 
laws. Josephson's law also holds, but at a rather sur- 
prising dimensionality which is the fractal one d ss 4.68 
[l9| . At the same time, the processes resulting in (3 ~ 1.9 
differ qualitatively from those resulting in (3 = 2, similar 
to thermodynamics where spatial fluctuations of the or- 
der parameter, neglected in mean field theories, result in 
j3 s» 0.33, therefore suggesting other universality classes 
were not completely ruled out. These could be the wave- 
wave interactions, responsible for the strong turbulence 
in plasma |20jp. which are next in dynamical importance 
and have fewer degrees of freedom [21| ■ 

In this paper, we use numerical simulations to study 
the threshold scalings in a weakly unstable collisionless 
Vlasov-Poisson system. Depending on the sign of the 
Poisson equation this set of equations describes either 
a plasma system or a gravitating system. The satura- 
tion mechanism in a collisionless gravitating system is 
the same as for the bump-on-tail instability in plasmas 
and threshold corresponds to the condition 7^ = in 
both cases. We show in Section [n] that the eigenfre- 
quency contains only an imaginary part, and therefore 
is the simplest model to study the threshold. Section llTTI 
describes the results of computations of the critical expo- 
nents and demonstrates that the scaling laws describing 
saturation are the same for plasma and gravitation. Sec- 
tion llVl addresses the scaling transformations of the phase 
space and the scaling law, which appears as a result of 
this symmetry. The exponent which describes correla- 
tions are obtained in Section [V] where Fisher's equality 
7 = i/(2 — 77) is also proved. In Section IVTI we show that 
the criticality in the system is described by the Dirac 
propagator for a fermionic field. We obtain hyperscal- 
ing laws and calculate upper critical dimensionalities in 
Section ETTI 



II. BASIC EQUATIONS 

The eigenfrequencies and eigenvectors of oscillations in 
a Vlasov-Poisson system are given by dispersion relation 



e[w(k),k] = , 



(9) 



where e is the permittivity (dielectric permittivity in the 
plasma case). The boundary between stable and unstable 
cases is determined by the condition 



Im[e(w,k)] = , 



(10) 
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|22| . For the bump-on-tail instability condition (|10fl sim- 
plifies to, Jl = Im(ijj) = 0, and criticality is related 
to the zero of the imaginary part of the eigenfrequency. 
Therefore we can employ a model which does not contain 
the real part; i.e., Rc(w) = 0. The simplest is the one- 
dimensional self-gravitating Vlasov-Poisson model which 
is described by the equations 



dt dx dx dv 
<9 2 $ 
dx 



"Hi2 = / f(x,v,t)dv - ! = 



(11) 
(12) 



where is the distribution function, and $ is 

the gravitational potential. Boundary conditions are as- 
sumed to be periodic in the x-direction. 
For the eigenfunctions 



X= X m exp(ifc m ir) 



m— — oo 



where 



X= [f(x,v,t),^(x,t)] T , 
X m = [f m (v,t),$ m (t)} T , 



(13) 



(14) 



(15) 



are the spatial Fourier components, the superscript T 
stands for transpose, k m — 2irm/L is the wavevector, 
and L is the system length we find 



fm "h ikmV fn 



m—m'-\-m 

Of, 



fm' d,V 



= 0, 



(16) 



or, explicitly for the components to = {0,1,2} and for 
L = 2n, 



d 

fo + ig^{Pif-i - P-ifi) = 0, 

d ( 1 

fi + iv j i + i-^ I pi/o + ^W-i - P-^h 

h + i + Qp 2 /o + Pi/iJ = 0, 
where 

/oo 
f m (v,t)dv , 
-oo 

is the Fourier component of density, and f-x = 
For a Maxwellian distribution 



(17) 
0(18) 
(19) 



A*. 



fo(v) = — ^exp 

the dispersion relation is 

1+ l<7 2 j(m) dZ(0 
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= 0, 
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FIG. 1: Amplitudes of the first four harmonics, m = 1, 2, 3, 4. 
The saturated level A sat is shown by the horizontal segment. 



where Z is the plasma dispersion function [23j, ( m = 
z m /\/2, z m = u) m / (k m <r), and a is velocity dispersion. 
In Eq. fify. 



2f \ 1 1 



(23) 



is the critical (Jeans) velocity dispersion for the mode to, 
and po is the background density. For small \z\ <C 1 [i.e., 
<7 2 /<7j(m) <C 1] the dispersion relation reads 



or 



£ [f^m i k m ) 
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a 2 j(m) 
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(24) 



(25) 



which is remarkably simpler than in plasma case, where 
the bump-on-tail instability and Landau damping ap- 
pear due to the wave-particle resonance at the phase ve- 
locity of the wave v p h = bJ pe /k. The frequency spectrum 
in the case of gravitation does not contain a real part, so 
the resonance occurs at v = 0; i.e., in the main body of 
the particle distribution. For all to, tr 2 (to) > cr 2 (to + 1), 
thus, if we write cr 2 (l) = <J 2 r , the distance from the in- 
stability threshold is 



(26) 



analogously to e = (T — T cr )/T cr . Using and J2HJ) 
the dispersion relation can be rewritten for uj\ as 



12 

7T 



(27) 



Time is measured in the units of the free-fall time 
tdyn = WGpo)^ 1 - Since AttG = 1 and the density 
po = m p n = 1, where m p is the particle mass, and n 



C=C„ 



is the concentration, tdyn 
here. 



2y / 7r in the units assumed 
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FIG. 2: The saturated amplitude A sa t vs. Aq for 9 — —0.05. 
Circles are calculated values, the curve is a spline approxima- 
tion. 

III. RESULTS 

Dispersion relation ll'ol) shows that there are no un- 
stable modes above the threshold a 2 r , a 2 > a 2 r , and 
therefore the system remains invariant with respect to 
translations x' — > x + r, where r is any number. Be- 
low the threshold the mode m = 1 becomes unstable, 
and therefore the continuous symmetry breaks and re- 
duces to a lower discrete one with respect to translations 
x' — > x + L. Therefore a 2 cr can be considered as the 
critical point of a second-order phase transition, and the 
amplitude of the mode m — 1 as the order parameter - 
following the definition of Landau . 

A. Order parameter scaling 

Equations and 112|) with initial distribution 

f(x, v, 0) = /r»[l + A cosihx)}, (28) 

were integrated numerically using the Cheng-Knorr 
method |2 1| . The amplitudes A m (t) = \p m (t)\ for m — 
1, 2, 3, 4 are shown in Figure^ The perturbation m = 1 
grows exponentially with the growth rate predicted by 
dispersion relation (|25p. Then, the growth saturates at 
some moment t = t sat at the amplitude A sat = Ai(t sa t)- 

Figure ^ a l so shows the (exponential) growth of per- 
turbations with m > 1 while i|25[l predicts exponential 
damping for these modes. This growth occurs because 
of nonlinear coupling between modes since the term pi/i 
dominates over the term P2/0 in equation (|19|l initially, 
when fi <C fi and therefore one has 72 = 2oji for the 
growth rate 72. 

Figure [21 shows that A sat is independent on Aq for 
small Aq, but there exists some threshold value of the 
initial perturbation Aq, when becomes dependent on Aq. 
This threshold amplitude A t hr corresponds to the trap- 
ping frequency uj = V 1 A t hr ~ At W(, ~ Wi the pro- 
cesses due to trapping become as important as of the 
resonance between wave and particles responsible for the 
linear Landau damping (or growth) in collisionless media. 
Therefore to rule out the influence of trapping processes 



FIG. 3: Amplitude A 3a t as a function of 6. Circles represent 
calculated data, the dashed line is the power law best fit. 

on linear growth the amplitude Aq must be small to pro- 
vide 

u b <wi- (29) 

The distribution function f(x,v,t) is plotted in Fig- 
ure 01 as a surface at the moment t = t sat . Note that the 
distribution function f(x,v,t) becomes flat in the part of 
the x — v domain 

v 2 

— + A sat cos(fcix) < A sat (30) 

separatrix, as predicted for the bump-on-tail instability 
[25ll2r1 |. Outside this area the Fourier component /i(w, t) 
remains modulated by the background Maxwellian dis- 
tribution as assumed at t = 0. and the components with 
to > 1 remain negligible [19J, so the dynamically im- 
portant area lies at v < \v sep \, where v sep = ±2\/A sat . 
The width of the dynamically important area must be 
small compared to a at maximum amplitude (i.e., A sat , 
Vsep c), otherwise the background distribution will be 
altered by evolution. 

Assuming the above two criteria, A sa t is calculated as 
a function of 8 and plotted in Figure |3 From Figure |3| 
we see that this dependence can be approximated by the 
power law 

A sat « ("<, (31) 

while (3 = 1.9950 ± 0.0034 for 6 < 0. Rewritten in terms 
of the bounce frequency lu and the linear growth rate 
7l = Im(u;i) the power law (|31|) becomes 

uj b = c-fP , (32) 

and the coefficient c = 3.22 ± 0.01. These values are 
almost identical to (3 = 2, and to the coefficient 3.2 in 
relation as predicted and calculated for the bump- 
on- tail instability USEE!. 

B. Response scaling 

Subjecting the system to an external pump of the form 
F(x) = F m cos,(k m x-\-(p) allows one to calculate the other 



FIG. 4: Distribution function (surface) and its isocontours 
(lines in the x — v plane) in configuration space at the moment 
of saturation t = t 3a t- 



two critical exponents, 7 and <5, which describe the re- 
sponse properties. The index 7 describes the divergence 
of the susceptibility, which can be written as 



dA sat {8) 



BFi 



(33) 



Fi ->0 



for m = 1. The results are shown by triangles (8 < 0) 
and circles (8 > 0) in Figure 

Computation of \ at some 8 requires at least 5 values of 
A sat corresponding to the given F\ . At the same time, F\ 
must small enough to avoid the effects of A sat depending 
nonlinearly on F\ . Again it requires extensive calculation 
of all quantities to high accuracy. In both cases x(#) is 
approximated by 



(34) 



and 7_ = 1.028 ± 0.025 for 8 < 0, 7+ = 1.033 ± 0.016 
for 8 > giving 7_ « 7+ = 7 « 1. These exp onents are 
very close to the corresponding results for [lJj , because x 
is the only linear coefficient, and this is common to both 
wave-particle and wave- wave interactions. 

The exponent 7 is the same as for the mean-field ther- 
modynamic models but, opposite to thermodynamics, 
the response is stronger at 8 < than at 8 > 0, as Fig- 
ure [S] shows. The susceptibilities are 



X- ~ 



(35) 



This difference, as well as the appearance of scaling 
and i|32[l with (3 = 2 instead of /3 = 1/2 can be ex- 
plained if one takes into account the difference between 
the Landau- Ginzburg Hamiltonian 



7~Llg 



A 
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(36) 



[where <f> is the order parameter and /i 2 ~ (T — T cr )/T cr ], 
which describes the Ising universality class, and the equa- 
tion 



y = y 
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FIG. 5: Plot of X vs. |0|. Circles (9 > 0) and triangles (9 < 0) 
represent experimental data, the dashed lines are the power- 
law best fits. 



which describes the amplitude of a weakly-unstable per- 
turbation in a one-species Vlasov-Poisson system [lfj| . 
According to (|37l) the maximum amplitude y sa t at y = 
scales with j L as y sat = t£. 

On the assumption that at 7l < the system responds 
linearly to an external pump dF, one can obtain the re- 
sponse dy to dF as 



3?; 2 

lL dy - -^fdy + 8F = 0. 



(38) 



However, equation l|37|l is not valid at 7^ < since it pre- 
dicts unlimited growth instead of damping in this case. 
At small initial perturbation yo the correct evolution is 
given by the linear equation y = ^Ly, and the suscepti- 
bility x = dy/dF is 



X+ = 1l ' 



and at 7^ > 



X- 



27Z 1 - 



(39) 



(40) 



At the critical point 8 = (or 7^ = 0) the response is 
described by another critical exponent S 



A sat oc F, 



1/6 



(41) 



The results of simulation are plotted in Figure El giving 
S = 1.503 ±0.005. This exponent cannot be obtained by 
the previous simple assumption from l|37|l because of its 
singularity at 7l = 0. 



IV. SCALING LAWS AND SYMMETRIES OF 
THE MODEL 

The remarkable property of the critical exponents 7, /?, 
and S is that they satisfy the Widom equality Q 14j with 
high accuracy. In thermodynamics the Widom equality 
is a consequence of the scaling of the Gibbs free energy 
under the transformation 



<8(\ a 'e,\ aB B) = \e(e,B) 



(42) 
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FIG. 6: A sa t as a function of F\. Circles represent experi- 
mental data, the dashed lines are the power-law best fit. 



from which it can be derived straightforwardly Q. The 
functions which comply the condition i|42|l are called gen- 
eralized homogeneous functions, and the condition itself 
is termed a homogeneity condition. 

The nature of this scaling for the marginally-stable 
Vlasov-Poisson system is clear from Figure 0] where the 
distribution function f(x,v,t) is plotted at the moment 
t = t S at- The remarkable property of the critical dynam- 
ics is the topological equivalence of the phase portraits 
for different 9: at the moments of saturation t\ and t 2 
corresponding to 9\ and 9 2 we can write 

f(x,X a "v,6 1 ,ti) = Xf(x,v,6 2 ,t 2 ) , (43) 



or 



f m (X a *v,6 1 ,ti) = Wm(v,6 2 ,t 2 ) , (44) 

for the Fourier component. Transformation between ti, 
6»i and t 2 , 9 2 can be written as t' — > X a *t, 9' — > X ae 9, so 



f m (X at t, X a "v,X a °9) = Xf m (t,v,9). 



(45) 



A weak external pump in the form F(x) = F% cos(kix+ 
ip) creates a similar topology in the phase space because 
of the same mechanism of the saturation and adds an ad- 
ditional variable to the distribution function. The trans- 
formation can be written as F[ — > A° F i F\ . Finally for 
fm = fm(ti v, 9, Fi) we can write the homogeneity condi- 
tion as 

f m (X a H, X a "v, X a ° 9, X aFl F x ) = Xf m (t, v, 6, F{). (46) 

The critical exponents (3, 7, and 5 can be expressed 
via the scaling exponents a v , ag, and ap 1 from which the 
Widom equality for the Vlasov-Poisson system can be 
proved directly (Appendix [SJ. They also provide a deep 
insight into symmetry properties of the system. Accord- 
ing to expression (|A4I) 



1 



P = 



rescaling the parameter of 9 (or growth rate) also rescales 
the distribution function f(t,v,x) in the ^-direction. 



This situation differs significantly from thermodynamics 
where 



1 - a B 



(47) 



This expression rescales the normalized distance from the 
critical point with external field B. 

Substituting A sat according to the power law (|31[1 for 
the order parameter to v^ ep = 4A sat oc (-0) 13 one can 
obtain (assuming a v = 1) 



(48) 



from which the scaling exponent is ag — 1 for (3 — 2 
(ag = 4 for — 1/2). Remarkably, the two different pro- 
cesses - the linear growth of an unstable perturbation due 
to the resonant wave-particle interaction and the subse- 
quent nonlinear saturation of this process due to particle 
trapping are interrelated. 

While there is no thermodynamic equilibrium in the 
collisionless system considered here, one can define the 
quantity which describes the response of the system to 
external thermal perturbation, just as the specific heat 
capacity describes the response of a thermodynamic sys- 
tem to heat transfer, C = 5Q/dT. For the case, consid- 
ered here 



C 



SQ 
d9 



dV_ 



(49) 



where V is the potential energy of the system. To cal- 
culate the specific heat capacity, V sa t corresponding to 
A sat is used. The critical exponent a can be calculated 
straightforwardly from H49fl and (|31fl . Because perturba- 
tions m > 1 are negligible for \9\ < 1, V sat oc A sat $ sat , 
where $ SQt = -A sat \ i.e., V sat oc - (- 9) 2/3 , < 0, and 
the heat capacity is given by 



where 



Coc-(-0) 



a = -(2/3-1). 



(50) 



(51) 



The scaling law (|5H can be proven using the homogene- 
ity condition (|46|l ( Appendix lA"|l . Remarkably, the criti- 
cal exponent a does not depend on the sign of Poisson's 
equation, and result is the same for the plasma case. 

Unlike thermodynamics where the relation between ex- 
ponents (3, 7, and a is given by Rushbrooke's equality, 
a + 2(3 + 7 = 2, the scaling law (|51|l does not contain 
the critical exponent 7. Nevertheless, the set of critical 
exponents a = -2.990 ± 0.006, = 1.995 ± 0.003, and 
7 = 1.031 ±0.021 satisfy Rushbrooke's equality with high 
accuracy. 



V. CORRELATION EXPONENTS 

The correlation function of fluctuations for the field 

dx 



E 



(52) 
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can be found from the fluctuation-dissipation theorem 
HI as 



(E 2 



T krx[e(w, k)] 
2txuj e 2 



(53) 



|2flj . where the permittivity e is given by l)24[l. Relation 
(|53[l can be integrated using the Kramers-Kronig disper- 
sion relations, and in the static limit — > (1531) becomes 



47TCT 2 

m p k B 



1 - 



1 



e(0,fc m ) 



(54) 



where k B is the Boltzmann constant and m p is the par- 
ticle mass. This equation can be rewritten as 



where 



47r a 2 



o 2 (rn) 



(55) 



aj(m) 



(56) 



The susceptibility \ can be written in terms of (E 2 



as 



X = (E 2 ) kl cx 0-\ 



(57) 



7 = 1. The combination of a and wi gives the character- 
istic length for the system from the dispersion relation 

(E3 



£ = 2tt- 



or, in terms of 9, 



£cx 



(58) 



(59) 



as (9 — > 0. Therefore the critical exponent that charac- 
terizes the correlation length is 1/ = 1. The correlation 
function [E 2 ) kl can be rewritten in terms of fcj = as 



(E 2 ) kl cx fc ; 



2- v 



(60) 



where 7/ is another critical exponent which characterizes 
the correlation function. On the other hand, using l|59l) 
one can rewrite this expression as 



(E 2 ) kl octf-^ 2 -") , 
and, taking into account (|57ll . 

from which finally we obtain the equality 
7 = v(2-rj). 



(61) 



(62) 



(63) 



The last equality is known as Fisher's equality and gives 
the last critical exponent, rj = 1. 



VI. RELATION WITH OTHER UNIVERSALITY 
CLASSES 

The correlation function l|55|) looks rather counter- 
intuitive, since at 9 rn > (damping waves), one has 



(E 2 



< 0, and the noise is imaginary. Nevertheless, 



this unusual situation has an analog - for particle-particle 
annihilation reactions of the type Y + Y — ► (corre- 
sponds to equation dn/dt — —an 2 , a > 0), Y — > 
(dn/dt — —an) the correlation function is also negative 
because of anticorrelation of particles |30j. In the case 
9 m > the amplitude A m — > as t — > 0. It is also shown 
that the criticality due to these annihilation processes 
belongs to a certain universality class which is different 
from the Ising universality class [30l l3l| and therefore is 
not described by the Landau- Ginzburg Hamiltonian |J5SJ. 

Another unusual quantity is the correlation length £ 
and the wavevector k^ = whose use allows us to 

establish the validity of Fisher's equality for the colli- 
sionless system, studied here. It is not related to the size 
of the system L but to the fluctuations in the system 
which determine an average path of correlated motion of 
particle in presence of these fluctuations. As the system 
approaches the threshold, fluctuations become correlated 
since the characteristic time of correlations ~ 9~ 1 
diverges as 9 — ► 0. This behavior is analogous to ther- 
modynamic systems where the correlation length is the 
only relevant scale near the critical point as e — > 0. 

To demonstrate that the criticality in the Vlasov- 
Poisson system belongs to a different class, let us compare 
the critical exponents corresponding to the Jeans insta- 
bility in a self-gravitating hydrodynamical system |27| . 
using the same approach. The dispersion relation for 
this system is 



(64) 



(65) 



where c 2 m = AnGp/k^ is the critical velocity of sound, 
corresponding to to 2 m = 0. As for kinetic case if c 2 > 
c\ = c 2 r , there are no unstable modes, and the correlation 
length is 
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2nc s 



-1/2 



(66) 



where 9f = (c 2 — c 2 r )/c 2 r is the reduced sound velocity in 
a fluid. Here we have the mean-field exponent Vf = 1/2. 

Assuming m = 1 and dividing the both sides of dis- 
persion relation l|65|l on c 2 one can obtain the correlation 
function as 



G^\h h ,9 f ) = 



-Of 



(67) 



This is the prop agator of Euclidean theory or of the scalar 
boson field 13] from which the Landau mean-field theory 
follows automatically. 
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On the other hand dispersion relation l|25[l for the col- 
lisionless case gives 



G (2) (fc ? ,#) = t 



2 k c 

7T fci 



(68) 



For collisionless systems the propagator thus corresponds 
to the vector fermionic field and describes a different 
class of critical phenomena. In the language of quantum 
field theory the parameters 8f and 8 are bare masses. 
Since 



G (2) (fc,0) oc 



1 



k 2- v ' 



(69) 



from l|H7| and (ffi%|l one can obtain i] = for the case of 
hydrodynamics and r\ — 1 for collisionless system. 



VII. HYPERSCALING LAWS 

The approach assumed in the previous section allows 
us to establish the hyperscaling law for the Vlasov- 
Poisson system which involves the dimensionality d along 
with critical exponents like Josephson's law ©■ Using 
propagator l|68[) which is the potential energy, the specific 
heat capacity C in d-dimensional space at — > can be 
obtained as 



C 



d_ 



d d k i G (2) (k^8) 



which gives 



c oc e~ d ■ 

With relation l(7T|l becomes 



Coc 8 



-v(2-d) 



(70) 



(71) 



(72) 



Taking into account the scaling law l|50|l for the specific 
heat capacity C one can obtain the hyperscaling relation 
which interrelates the exponents a, v, and the dimen- 
sionality d 



a = v{2~ d) 



(73) 



The last equality reveals d = 2 as the upper critical 
dimensionality for the Vlasov-Poisson system since the 
heat capacity becomes divergent if d < 2, thus indicating 
the importance of fluctuations in the critical area. It also 
shows that the dimensionality corresponding to the crit- 
ical exponents a = — 3 and v = 1 is d = 5, fluctuations 
at 8 w arc insignificant, and therefore a — — 3, (3 = 2, 
7 = 1, v = 1, and r/ = 1 are the mean-field exponents. 

The use of the scalar field propagator l|67|l instead of 
gives 



1/(4 - d) , 



(74) 



and at a — the upper critical dimensionality is d c = 4 
which is the Landau mean-field theory case for the Ising 



universality class. However, relation l|74|l is not valid for 
the Vlasov-Poisson system because of its different prop- 
agator. On the contrary to relations l|73ll and (|74|l which 
are valid for specific propagators l|67l) and H68|) . Joseph- 
son's law is universal for all cases considered. With 
exponents v = 1 and vj = 1/2 it gives d c = 2 and d c = 4 
as the upper critical dimensionalities for the collisionless 
and hydrodynamic cases, respectively, and d = 5 for the 
exponents of the Vlasov-Poisson system calculated here. 
Without going into details here, we note that this uni- 
versality appears because the fundamental description is 
given by the same functional integrals in both cases. In 
particular for the free scalar bosonic field (no interac- 
tions) the partition function is 



Zq = / T>4> exp 



d d xH 



where 7io is the Landau-Ginzburg Hamiltonian Hlg <l3l)|) 
without quadratic term. In the fermionic case the La- 
grangian for a Dirac spinor field is used instead of Ho- 



VIII. CONCLUSIONS 

We have studied numerically and analytically a model 
Vlasov-Poisson system near the point of a marginal sta- 
bility. The most important finding is that the critical- 
ity of the Vlasov-Poisson model studied here belongs to 
a universality class described by the propagator corre- 
sponding to a fermionic vector field. This finding is in 
striking contrast with the previous critical phenomena 
studies concerning systems whose criticality belongs to 
universality classes corresponding to the scalar bosonic 
fields, like the Ising universality class. 

This fundamental discrepancy emerges from the 
qualitative difference between objects considered: the 
Landau-Ginzburg Hamiltonian l|36l) takes into account 
spatial variations of the order parameter via the local 
differential operator V, whereas the integro-diffcrcntial 
operator for the Vlasov-Poisson model acts on the dis- 
tribution function containing the additional dimension of 
velocity. 

We have calculated numerically the critical exponents 
which describe the critical state of the model and estab- 
lished analytically that these exponents and the dimen- 
sionality are interrelated by the scaling and hyperscaling 
laws like the Widom, Rushbrooke, and Josephson laws 
at the formal dimensionality d — 5. The upper critical 
dimensionality is d c = 2 and since d > d c the calculated 
exponents are the mean-field exponents, different from 
those which one might expect the Landau- Weiss set of 
critical exponents corresponding to the Ising mean-field 
model where d c — 4. This is related to the higher di- 
mensionality of the Vlasov-Poisson kinetic problem as- 
sociated with the velocity space and to the type of the 
criticality of the Vlasov-Poisson systems, which belongs 
to a universality class different from the Ising universality 
class. 
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The critical exponents we have found here are a = —3, 
f3 = 2, 7 = 1, 8 = 1.5, v = 1 and r\ = 1. The difference 
between this set and the set a « —2.814, (3 « 1.907, 
7 » 1, 5 w 1.544, u = 1 and 77 = 1 [lj| is because ^4 sat 
is about 50 times larger for the latter case, thus caus- 
ing wave-wave interactions to dominate, thereby yield- 
ing a different universality class. More important, the 
later exponents satisfy scaling laws at fractal dimension 
d fts 4.68 indicating reduced dimensionality because wave- 
wave interactions have fewer degrees of freedom than 
wave-particle ones |2l| . 
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APPENDIX A: RELATION BETWEEN THE 
SCALING AND CRITICAL EXPONENTS 

From the homogeneity condition 1|46|) 

f m (\ at t, X a "v, \ ae 9, X aFl F 1 ) (Al) 

= \f m (t,v,e,F 1 ) , 



Equations (|A4|) - (|A6p can be rewritten in matrix form 



for p m components by integration over v, one has 

A~°" p m {X at t, \ a °6, X aFl Fx) 
= Xp m (t,6,Fx). 



(A2) 



For any two A sat = pi{t sat ) and A' sat = pi(t' sat ) one can 
write 



as 



WA = X 



(A8) 



where A = [ag, a v , a Fl ] T , X = [1, — 1, — S] T , and the 
trix W is 




(A9) 



(A10) 



Using l|A6(l to eliminate a v , the system (|A8ll can be 

reduced to 



W = 

The determinant of W is 

detW = -/3 + (S/3-7 = 



flag - -a Fl = , 
a 



lag 



1 )a Fl = , 



(All) 
(A12) 



for which solution exists only if the Widom equality 7 = 
P{5 — 1) holds. Therefore ag and a Fl can be formally 
considered as the eigenvectors of W whose eigenvalue is 
A = 0. In particular 



1 



/3 + 7 



a Fl 



(A13) 



which indicates that rescaling of the distribution function 
under an external pump is equivalent to rescaling due to 
the field which appears for nonzero order parameter. 



X- a *A sat {X ae 0, X a ^Fx) = XA sat (6,Fx). (A3) 

Assuming A = (— 1/0) 1 / ", the critical exponent (3 can 
be rewritten in terms of the scaling exponents a v and ag 
as 



1 + a v 
ag 



In the similar way for 7 and 5 one can write 
-a v — 1 + a Fl 



7 



ag 



5 = 



l + a v 



(A4) 



(A5) 



(A6) 



and the Widom relation follows from (|A5() straightfor- 
wardly: 

_ — a v — 1 + a Fl _ 1 + a v a Fl 



a 



-13 + 138 
(3(8-1) . 



(A7) 



APPENDIX B: RUSHBROOKE'S LAW FOR 
VLASOV-POISSON SYSTEM 

The heat capacity can be formally defined as 



C= SQ = dV 
~ dO ~ d6' 



(Bl) 



where V is the potential energy of the system. To cal- 
culate the specific heat capacity, V sa t corresponding to 
A sat is used. 

Because perturbations m > 1 are negligible for \6\ <C 1, 
V sat oc A SQt $ sat , where $ sat = -A sat , and 



V sat oc A sat . 
From HA3fl one can obtain 



(B2) 



^X- 2a *A 2 sat (X a °6, X a ^Fx) = -^X 2 A 2 aat (9,Fx), (B3) 



or 



^X- 2a "- 2 A 2 sat (X a °e, X a ^Fx) = ^A 2 sat (6,Fx) . (B4) 
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Assuming A = 6 1 / a " and F\ = 0, Equation i|B4|) can Equation (|B5|) has the form of the power law, C(9,0) oc 
be rewritten as 8~ a , with 

^[fl( J - +8 )/"»4(-l l 0)] = ^A 2 sat (9,0), (B5) 

or a + 1 

a = _ 2 o 1 +i + 1 = _ 2(3 + 1 (Bg) 

AL(-i,o) = |aL(^o), 

ag of) 

(B6) 

or 

o i o (2 +2> The last relation corresponds to Rushbrooke's equality 
A^ t (-l,O)0^ 1 =(7(0,0)- (B7) a + 2/3 + 7 = 2 at 7=1- 
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